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Abstract: We report on initial results from a project to constrain the large-scale and turbulent magnetic 
fields of the Milky Way galaxy, which eventually will incorporate all of the relevant observational data. 
In this paper we fit popular large scale magnetic field models to WMAP3 polarization maps. We find that 
the polarization data can constrain certain model parameters but does not uniquely determine the best-fit 
parameters. We also find that the polarization data alone cannot distinguish between model symmetries, 
e.g., the existence of field reversals. We show how future UHECR data can break this degeneracy. 
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Introduction 

A variety of observational methods has been used 
to constrain the magnetic field of the Milky Way, 
e.g., synchrotron emission, Faraday rotation. Zee- 
man splitting and polarization of starlight. In 
a project currently underway also incuding M. 
Haverkorn and A. Lazarian we make use of all rel- 
evant data available to study the Galactic magnetic 
field (GMF), at large and small scales. The rea- 
son for combining different types of data goes be- 
yond having a larger number of data points when 
performing parameter estimation. Different types 
of data constrain different aspects of the magnetic 
field. For instance, synchrotron emission probes 
the integrated magnetic component perpendicu- 
lar to the line-of-sight while Faraday rotation and 
starlight polarization probe the parallel component. 
Zeeman splitting allow in situ measurements of the 
magnetic field. 

In this paper we present first results, using po- 
larization data from the Wilkinson Microwave 
Anisotropy Probe (WMAP) 3rd year data release 
[3]. A previous analysis of fitting the WMAP3 data 
to a large scale magnetic field model was done in 
[3]. There the WMAP team fitted the polarization 
angle of synchrotron emission as predicted by an 
axisymmetric magnetic field to the measured po- 
larization angle (j—j tan^' (f/ / Q)), after masking 



out 25.7% of the sky (essentially the Galactic disk 
and the North Galactic Spur). 

Galactic magnetic field models 

When describing models of the regular part of 
Galactic magnetic fields it is natural to use cylin- 
drical coordinates (r, 9,z). Typically, the disk field 
is then parametrized by 

B,.^B{r,Q)cosp, Bq = B{r,d)smp, (1) 

where p is the pitch angle. There has been some 
confusion in notation in the literature, and we will 
use the following conventions. A field that obeys 
B{r,Q) — B{r,Q + a) for any a is called axisym- 
metric, if B{r,Q) = — Z?(r, 9 + 7t) it is called bisym- 
metric, and if B{r,Q) = B{r,Q + k) we propose to 
call it disymmetric. Another distinction that can 
be made is a model's symmetry properties under 
z ~z, i.e., reflection at the disk plane. We 
call a field symmetric with respect to the Galactic 
plane if B(r, 9,— z) =B(r, 9,z), and antisymmetric 
if B(r,9,— z) = — B(r,9,z). This notation agrees 
with, e.g., [6, 2], but conflicts with, e.g., [5]. 

A standard way to parametrize a bisymmetric spi- 
ral field [1] is 

fi(r,9) =/?(r)cos(9-pin-^), (2) 
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where 9 is defined to increase clockwise around the 
Galactic center, and 6 = points to the Sun. The 
pitch angle p is positive if the clockwise tangent to 
the spiral is outside a circle of radius r centered on 
the Galactic center, and (3 = 1 / tanp. At the point 
(ro, = 0°), the field reaches the first maximum in 
the direction I — 180° outside the solar circle. 

We choose b{r) = bo for r < rc, and b{r) = bo^ 
for r > Tc, and define the field for r < 20kpc. 

Using WMAP3 polarization data 

Relativistic electrons accelerated by a magnetic 
field will radiate synchrotron radiation [4]. For 
a power law distribution of cosmic ray electrons, 
n(E)dE ~ E~^dE, the synchrotron emissivity is 

1+1 i-i 

jv^fi^v^. (3) 

The emission has a large degree of Unear polariza- 
tion. 

WMAP is an all-sky survey measuring the Stokes 
I,Q,U parameters in five frequency bands in the 
23-94 GHz range. The Stokes / parameter is pro- 
portional to the total (polarized and unpolarized) 
emission, and the Q and U parameters describe 
the polarized part of the emission. In the K band 
(23 GHz) the measured polarized radiation is dom- 
inantly Galactic synchrotron emission. 

We obtain the best fit by minimizing for Q and 
U. The resolution in the K band is ~ 1°, but to bet- 
ter probe the large scale regular field we smooth 
the data to ~ 4°. For each 4° pixel we calculate 
the variance of the 1° pixels within a radius p from 
the center position of the larger pixel. This way 
we obtain Og and Ou for the 4° maps. Regions 
with large turbulence and irregular structures have 
a larger variance and thus are deweighted in the 
fit. We thus avoid masking out the Galactic disk, 
which is essential to be able to probe the magnetic 
structure in the plane of the Galaxy. For a given 
Galactic magnetic field and electron distribution 
(using the Hammurabi code [8]) we calculate the 
polarized synchrotron emission Q and U. We then 
measure the goodness-of-fit by 

pixels /^nwdel ndata\2 

il=L ^^' ) , (4) 

(=1 "Q,( 



and similarly for U, with Xq^ =y^-\-y^ the quan- 
tity that is minimized. 

In figure 2 a smoothed map of the WMAP3 Q pa- 
rameter is shown, along with map of Oq obtained 
using p = 2°. 

Results 

We perform a parameter search to minimize Xq[/ 
for the parameters bo, and zo for a bisymmet- 
ric spiral field that is symmetric in z- Following 
Stanev [5], we put p = —10° and ro = 10.55kpc. 
From figure 3 it is evident that the data strongly ex- 
clude large regions of parameter space. However, 
no unique best-fit parameters are found, as a num- 
ber of combinations of the search parameters yield 
a x^ very close to the minimum. 

To investigate the effect of model symmetries, we 
perform the same analysis on a bisymmetric field 
that is antisymmetric in z, and two disymmetric 
field models (symmetric/antisymmetric in z). We 
obtain the disymmetric model from equation 2 by 
squaring the cosine. We find that the minimized 
is almost identical for all four models. 

To examine further the hmitations of using polar- 
ization maps to constrain GMF models we pro- 
duce mock Q and U maps with Hammurabi us- 
ing a z- symmetric bisymmetric field model with 
the parameter values bo = 3.5;uG, Vc = lOkpc, 
zo = l.lSkpc, p = -10° and ro = 10.55kpc. Us- 
ing these maps as "data", and with uniform a, we 
do a parameter search for the same bisymmetric 
field used to obtain the mock data. The result is 
shown in figure 4. We note that the correct best- 
fit values are recovered, but that the insensitivity to 
the model parameters remain. If a z-antisymmetric 
disymmetric model is fitted to the above mock data 
we find the surfaces to be visually indistinguish- 
able to those in figure 4, and we can conclude that 
even with ideal data, polarization maps alone can- 
not distinguish between the various model symme- 
tries. 

UHECR multiplets 

A powerful way to break this degeneracy between 
model synmietries would be to use future UHECR 
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Figure 1 : A hypothetical UHECR multiplet source 
along with the UHECR locations as predicted from 
bisymmetric and disymmetric spiral GMF models. 
Subscripts refer to symmetry/antisymmetry under 
reflection in z- UHECR energies selected in the 
lO^^-^OEeV range, with the size of the circles pro- 
portional to the event energies. 



multiplets. UHECR deflection is unique among 
observables in being independent of electron den- 
sity and very weakly dependent on random fields 
because the correlation length of the random fields 
are much smaller than the Larmor radius [7]. As 
the deflection angle of a UHECR is inversely pro- 
portional to its energy and proportional to the 
transverse magnetic field, the arrival directions of 
an UHECR multiplet will roughly form a string on 
the sky; the highest energy cosmic ray closest to 
the direction of the source, and the lower energy 
cosmic rays further removed according to their en- 
ergy. From their energies and angular position on 
the sky their common source can be estimated. A 
hypothetical multiplet placed in the southern hemi- 
sphere is shown in figure 1, together with the ar- 
rival directions of the multiplet as predicted for the 
best fit spiral field parameters with the four differ- 
ent GMF symmetries. It is clear that if a number 
of such multiplets would be discovered in future 
experiments they would be very useful in breaking 
the degeneracy of the model symmetries. 



Summary 

We have reported on a project that will eventually 
make use of all relevant data to constrain models of 
the Galactic magnetic field. In this paper we out- 
line our analysis using the WMAP3 polarization 
data. We improve on the previous analysis [3] of 
the same data by including the Galactic plane, and 
allow for weighting individual pixels based on the 
variance in the data. 

We find that polarization maps alone can exclude 
large regions of parameter space for GMF mod- 
els common in the literature, but fail to give a 
unique set of best-fit parameters. Furthermore, we 
find the data cannot distinguish between funda- 
mental model symmetries, such as the existence of 
field reversals between spiral arms. We comment 
on the usefulness of using UHECR multiplets to 
break this degeneracy, if such multiplets were to 
be found. 

In future work random fields will be introduced 
and the synchrotron component of /, in addition to 
Q and U will be used in the fit, as well as Faraday 
Rotation Measures and starlight polarization data. 
More details will be given in a forthcoming paper 
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Figure 4: surface plots for bisymmetric field model fitted to mock data from the same bisymmetric GMF 
model, values only shown for relative comparison. 



